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Introduzione 


Questo note si propone di portare all'attenzione degli speleologi 
che si apprestano a fare misurazioni di tracciamenti aerei alcune 
semplici formule per fare una primaria analisi quantitativa dei 
risultati. 


Aria 


Gia' oltre un secolo fa' H. Bock face una modellizzazione fisica 
delle differenti circolazione dell'aria nelle grotte, distiguendo 
quelle che chiamiamo "cavita' barometriche" e "cavita' a tubo di 
vento". Nella figura si vedono chiaramente distinte la 
circolazione invernale e quella estiva, ad indicare una grotta con 
circolazione a "tubo di vento", ovvero una cavita' in cui un 
ingresso inferiore e' collegato ad uno superiore. 

Il fatto che l'osservazione dei flussi dell'aria nelle grotte possa 
dare informazioni utili alla ricerca ed esplorazione speleologica e' 
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Matbematifch-phyfikalifche Unterfuchung der 
Eisbòblen und Windròbren. 


Von Ingenieur Hermann Bodk. 


ormai assodato da decenni, e spesso si tende ad insistere di piu' 
su una strettoia ventosa che su una dove non c'e' (flusso di) aria 
[Badino, Maistre]. 

Dato che un flusso di aria fra due ingressi e' indicazione di una 
connessione fra essi, anche se poi, non esiste (ancora) un 
passaggio praticabile da alcuni anni gli speelologi hanno iniziato 
a fare misurazioni che permettessero loro di identificare 
connesioni e quindi possibilita' esplorative [Miragoli]. 


Cosa si puo' misurare agli ingressi ? Prima di tutto le portate. 
Tanta aria deve uscire quanta ne entra. Se da un ingresso alto 
esce meno (massa di) aria di quanta ne entra dall'ingresso 
basso vuol dire che parte dell'aria che entra esce da qualche 
altro ingresso. Pero' la misura della portata non e' affatto 
semplice, perche' le sezioni degli ingressi sono sovente irregolari 
e la velocita' del flusso non e' uniforme su tutta la sezione. Una 
soluzione per una stima accettabile e' suddividere idealmente la 
sezione in rettangoli/triangoli, e misurare la velocita' al centro dei 
rettangoli/triangoli. Infine sommare i prodotti area X velocita'. Il 
problema principale e' la misura della velocita' (v. figura sotto). 
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Area(A)=150x30 cm? 
Area(B)=90x30 cm ? 


Velocita'(A) = 20 cm/s 
Velocita'(B) = 10 cm/s 


Portata = 0.45 x 0.2 
+0.27 x 0.1 
= 0.12 mÌ/s 


particolare la temperatura. L'analisi delle variazioni giornaliere e/ 
o stagionali di questa permette di fare correlazioni fra i vari 
ingressi ed ipotizzare connessioni. 

La misura dei parametri chimici (anidride carbonica, anidride 
solforosa, ...) dovuti a cause naturali risulta rilevante solo in 
particolari situazioni, oppure per rilevare problemi ambientali di 
inquinamento. 

Diventa invece importatante quando si fanno tracciamenti aerei. 


Tracciamenti aerei 


| tracciamenti aerei consistono nell'immettere nell'aria (di un 
ingresso aspirante, oppure in un particolare punto all'interno 
della grotta) un composto volatile e di posizionare sensori con 
datalogger in altri punti di rilevazione. In genere l'immisione del 
tracciante ha durata molto breve, idealmente "istanatanea". 
Dall'analisi delle curve di risposta si possono trarre informazioni 
ai fini esplorativi. La sostanza immessa deve essere volatile, 
mesclarsi facilmente con l'aria, non esser adsorbita all'interno 
della grotta, e non risultare nociva per la salute umana, ne' 
dannosa per l'ambiente. 

Questo approccio ha avuto sviluppo nell'ultimo decennio. Sono 
state provate diverse sostanze traccianti, terpeni [Cella 2015], 


Tipica situazione di un 
tracciamento aereo. La 
porzione coperta dal quadrato 
nero rappresenza una parte 
di grotta non conosciuta di cui 


si vuole verificare l'esistenza. 


anidride carbonica [Cella 2019], e VOC (Volatile Organic 
Compound). 

Una differenza fra i tracciamenti aerei e quelli idrologici e' che 
l'aria scorre nelle grotte sia verso il "basso" (circolazione estiva) 
che verso l'alto" (circolazione invernale), Quindi e' possibile 
scegliere il verso del tracciamento adeguandosi alla stagione. 
D'altra parte le stagioni intermedie, cioe' i periodi di inversione 
non vanno bene per i tracciamenti aerei. 


Che cosa si misura? 

La quantita' (grammi) di tracciante immesso, e l'orario di 
immissione. Si misura poi la curva di restituzione, solitamente 
espressa in ppm (parti per milione) o ppb (parti per miliardo) in 
funzione del tempo. | valori della curva di restituzione sono 
solitamente campionati (cioe' misurati) ad intervalli regolari, 
tipicamente poche decine di secondi. 

La figura sotto riporta un esempio di curva di restittuzione. Sono 
riportati i valori (ppm) in funzione del tempo. Da notare che il 
valore di base non e' sullo zero della scala verticale. 


Curva di restituzione. A. Vernassa 2023. 


Analisi qualitativa 


Come prima cosa si puo' fare una analisi qualitativa. 

Prima di tutto la presenza o meno di un collegamento (aereo): in 
un tracciamento aereo si immette il tracciante in un sol punto ma 
si mettono i rilevatori in piu' punti dove ci si aspetta di rilevare il 
tracciante, oppure si ci si aspetta di non-rilevarlo. 

La forma della curva di restituzione e la sua regolarita' danno 
delle informazioni sul tipo di percorso seguito dall'aria. Una 
forma semplice come quella alla pagina precedente e' indice di 
un solo percorso fra punto di immissione e punto di rilevazione. 
Una forma con due picchi indica la presenza di due percorsi 
alternativi, entrambi seguiti dall'aria. Una forma molto irregolare 
e distribuita su un longo tempo viene interpretata come 
ostruzioni tipo frane longo il percorso dell'aria. 
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Esempi di curve di restittuzione. 
G.D. Cella - Bossea 2021 


Analisi quantitativa 


L'analisi quantitativa basata sui momenti della curva di 
restituzione permette di ottenere stime della portata e del volume 
della grotta. Questo approccio e' stato proposto in campo 
medicale e risale a oltre 50 anni fa [Zierler]. Per una 
formulazione rigorosa e per la dimostrazione rimando all'articolo 
di Lightfoot et al. In particolare l'analisi si basa su tre ipotesi 

e non c'e' diffusione del tracciante attraverso le pareti 
della cavita'; 

e non c'e' assorbimento del tracciante; 

* l'intervallo di immissione e' molto corto; 

* la velocita' di diffusione e' trascurabile rispetto al moto 
convettivo all'ingresso e all'uscita. 
Le prime due ipotesi sono abbastanza plausibili, anche se 
nell'attraversamento di grandi ambienti il tracciante potrebbe 
esser "disperso" in sacche d'aria stagnante. La terza ipotesi e' 
accettabile. E' possibile tener in conto l'intervallo di immissione, 
ma in pratica l'effetto e' piccolo (1%) e solitamente trascurabile 
nell'analisi. Le velocita' del moto convettivo sono dell'ordine di 
alcuni cm/s o dm/s, mentre la diffusivita' delle sostanze usate e' 
di circa 10 mm?/s. La velocita' della diffusione e' proporzionale al 
gradiente della concentrazione. Quindi, ci deve esser una 
discreta corrente d'aria nei punti di immissione e rilevazione. 


Il momento di ordine n della curva di restituzione e' definito come 
M© = fr(t) t" dt 

dove il tempo e' misurato dall'istante di immissione, e r(t) 

rappresenta il valore in ppm (o ppb) sulla curva di restituzione, 

cui e' stato sottratto il valore di base 


Per esempio per calcolare il momemto di ordine uno, M&), 

usando i dati sulla curva a pagina 5, devo fare la somma di 
(R(t) - Ro) (t-t) Dt 

dove R(t) e' il valore del punto sulla curva, R; il valore base (che 

si evince dalla parte iniziale della curva a sinistra), t; il tempo di 

immissione, a Dt la differenza di tempo fra due valori campionati 

consecutivi. 

Da notare che se i valori sono in parti per milione bisogna 

dividere per 105, se in ppb bisogna dividere per 10°. 

Alla fine il risultato e' espresso in [parti-S7]. 


Si puo' dimostrare che la portata, Q, al punto di rilevazione e' 
data dal rapporto fra la massa di tracciante immessa e il 
momento di ordine zero, e il volume, V, della grotta e' dato dal 
rapporto fra il momento di ordine uno e quello di ordine zero. Per 
la dimostrazione rimando ai due lavori citati in bibliografia. 


Tenendo conto delle unita' di misura risulta 
Q = (Valm?] / Ma[g] ) m[g] / M9so[s] 


Va volume molare (RT/P) 
Ma massa molare tracciante 


VE Q[m*/s] MA)s2[S7] / M9s2[S] 


Da notare che la portata dipende dalle condizioni del gas: 
R=8.31 J/g-mol-K, la temperatura (in gradi K) e la pressione nel 
punto di rilevamento. 


A titolo di esempio queste formule sono state applicate al 
tracciamento effettuato da A. Vernassa nella grotta Buranco di 
Bardineto (la cui curva di restituzione e' riportata a p. 5). La 
figura sotto e' l'ingresso dove e' stato posto il rilevatore. La 
sezione e' di circa 1 m2. Sono stati immessi circa 500 g di 
tracciante propano intorno alle 11:50 (10 minuti prima delle 12). 
La tabella sotto riporta alcuni dati della curva di restituzione. | 
valori sono campionati ad intervali di 30 s. Il valore base e' 29.7 


ppm. 


Buranco da Sabbia 


I termini nel calcolo di M® corrispondenti a questi due dati 
saranno 

148.6 x 30 / 1000000 parti-s 

180.8 x 30 / 1000000 parti-s° 
Quindi per calcolare M basta fare la somma dei valori nella 
quarta colonna e moltiplicare per 30/1000000. 


Similmente quelli nel calcolo di Ml) 

33861 x 148.6 x 30 / 1000000 parti-s 

33891 x 180.8 x 30 / 1000000 parti-s 
Pertanto bisogna fare i prodotti della terza e quarta colonna, 
sommarli, e moltiplicare il risultato per 30/1000000. 


Come volume molare supponiamo 0.023 m3/mol, e la massa 
molare del propano, C,Hj, circa 44 g/mol. Sommando tutti i 
valori risultano all'incirca 

MO) = 5.72 s 

MM= 3.64 105 s2 
da cui si trovano 

Q = 0.046 mì / s, cioe' una velocita di circa 5 cm/s 

V = 2800 m3 
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Considerazioni 


Come abbiamo detto la misura della portata e' difficile. Tuttavia 
e' possibile che il valore stimato differisca significativamente da 
quello osservato anche tenuto conto della grande incertezza di 
quest'ultimo. In questo caso si possono trarre interessanti 
conclusioni. 
Infatti noi rileviamo il flusso fra un punto di immissione ed uno di 
rilevazione, ma non e' detto che questo sia l'unico percorso 
dell'aria. Ci possono essere due ingressi connessi all'uscita. In 
questo caso il flusso osservato all'uscita e' la somma dei flussi 
provenienti dai due ingressi percio' 

Q(stimata) < Q(osservata) 


Quando invece c'e' una seconda uscita il flusso in ingresso si 

distribuisce sulle due uscite e pertanto solo parte della massa 

immessa arriva all'uscita col rilevatore, mentre la stima della 

portata considera tutto il tracciante immesso. Quindi 
Q(stimata) > Q(osservata) 


Se la portata stimata e' in accordo con quella osservata, non 
possiamo trarre alcuna indicazione riguardo altri ingressi o altre 
uscite. 
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Derivazione delle formule 


Si suppone di avere un flusso costante attraverso una cavita' di 
forma arbitraria, con un solo ingresso ed una sola uscita (figura 
p. 7). Si immette instantaneamente una massa m di tracciante 
all'ingresso. La distrubuzione del tracciante e' descritta dalla 
equazione di diffusione 


Op/ot = - (D - pv) + (D :dP DP) 


dove p e' la concentrazione locale del tracciante (massa per 
unita; di volume), 2p la diffusivita' del tracciante in aria, v la 
velocita' dell'aria, t il tempo, e D e' l'operatore di derivazione. Se 
il moto e' turbolento si considerano le quantita' mediate sul 
tempo e la diffusivita' effettiva per moto turbolente. 
Si assume che non c'e' diffusione attraverso le pareti della grotta 
Sf 

n:v=0, n-Dp=0, 
e che il flusso diffisivo del tracciante sia piccolo in confronto al 
flusso convettivo all'ingresso e all'uscita (superifici S, e S,) 

De (Nn: Dp)<<p(n-v). 
Il vettore n e' il vettore unitario normale alla superficie, uscente 
(figira p. 7). 
La concentrazione del tracciante all;ingresso S1 e' zero fino al 
tempo t=0, e dopo il tempo t=t,. In pratica l'impulso di tracciante 
dovrebbe avere una durata molto breve. 


| momenti rispetto al tempo della concentrazione sono 

Ir) = f p(r,t) t" dt 
dove l'integrale va da -c0 a +00. Moltiplicando per tn l'equazione di 
diffusione ed integrando sul tempo (limitandoci ai valori in cui la 
concentrazione in uscita e' non-zero, cioe; fra 0 e 00) ed 
integrando per parti 

fop/attr dt=-(D-vfpt"dt)+(D-2Dfpt" dt) 


=[pt"- nfpt-1dt 

Quindi si ha un insieme di equazioni 

-n105=-(D-vI®)+(D -DDIM) 
con il termine a sinistra zero per n=0. Integrando queste 
equazioni sul volume della grotta ed applicando il teorema della 
divergenza di Gauss si ha, per n=0, 

=-f(n:vI®)dS+f(n-DDIO) ds 

L'integrase su Sf non da' contributo per le condizioni al contorno 
e negli integrali sulle superfici di ingresso ed uscita si puo' 
trascurare il termine diffusivo. Quindi 

0=fv IO dS1-fv 10 dSs2 
e supponendo I(0) costante sulla sezione di ingresso e su quella 
di uscita si ha 

<v> Si 1,09 = <v> S, 109 
Il prodotto <v>S e' la portata Q. Quando questa e' costante 1,9 
= 1,09, cioe' 109 e' uguale sul piano S, e sul piano S,, e per 
l'arbitrarieta' del piano S,, possiamo dire che 10 e' costante in 
ogni punto della grotta. Il valore di 19 viene denotato M,, 
momento di ordine zero. Calcolando 1° sul piano S,, si ha 

1,9 = p(r,t) dt = f m/V è(t) dt = m At/V = m/Q 
dove l'integrale fra 0 e co da' il tempo perche' il tracciante passi 
nella grotta. In conclusione si ottiene 

Q=m/M 
che permette di stimare la portata conoscendo la massa di 
tracciante immessa e calcolando il momento di ordine zero 
all'uscita. 


integrando sul volume della grotta la relazione per gli altri valori 
di n, applicando il teorema della divergenza per l'integrale dei 
termini a destra e tenendo conto che il termine diffusivo e' 
trascurabile rispetto a quello convettivo sulle superfici di 
ingresso ed uscita, 

- n (1/7) f 10 dV = - (1/V) f v 10 ds 
e nell'ipotesi che I(n) sia costante sulle sezioni di ingresso e 
uscita 


14 


15 


- n (1/V) f 101) dV = (1/V) ( <vj> S, 10 - <v,> S, 1,10) = 
(Q/V) ( I, 7 1,0) ) 
Se la concentrazione del tracciante sulla superficie S, e' una 
funzione delta (cioe' immesso in un intervallo molto breve) 1, = 
0 e questa equazione, per n=1, permette di stimare il volume 
della grotta 

V=Q1,05)/M, 


Anche momenti di ordine superiore sono utili, in particolare il 
secondo momento centrale 

Hu, =Sp(t-t*)? dt 
dove t* e' il tempo medio di residenza del tracciante, t* = V/Q. La 
varianza 02 = |, / (M.)? descrive l'ampiezza della curva di 
restituzione. 


Pressione e temperatura 


La temperatura nel punto di rilevamento dovrebbe restare 
abbastanza stabile durante l'acquisizione, essendo la 
temperatura dell'aria in uscita dalla grotta e quindi poco 
influenzata sul breve periodo da variazioni giornaliere esterne. 
Una variazione di 1°C corrisponde ad un errore di 0.3%. A tal 
fine e' bene che il punto di rilevamento sia posto un poco 
all'interno dell'ingresso dove esce l'aria. 

La pressione invece risente delle variazioni di pressione esterna. 
Quando il tempo e' stabile si tratta di variazioni di pochi mbar, 
quindi un errore dello stesso ordine. In caso di tempo instabile la 
pressione puo' variare di decine di mbar, quindi non tenerne 
conto si traduce in un errore di qualche percento. 


Volendo tener conto delle variazioni di pressione e temperatura 
bisogna includere il volume molare nel calcolo del valore MO 
usato per la stima della portata. Quindi si deve usare 

M'0) = 1/R f r(t) P/T dt 
Per la stima del volume si usa MO). 


Interpolazione lineare dei valori campionati 


Nel calcolo dei momenti si moltiplica il valore della curva di 
restituzione in un dato istante per l'intervallo, D, di tempo di 
campionamento. Si potrebbe pensare di sommare per ciascun 
secondo i valori interpolati (lineramente) fra due istanti di 
campionamento. 
Ovvero 
n Leo 3 f(m) + i/D (f(n+1) - f(mM)]= 
=> [Df(n)+(D-1)/2 (f(n+1) - f(n))] 
=D f(n) 
poiche' il secondo termine e' zero giacche' f(n)=0 agli estremi 
della somma. 
Similmente 
In Zio II f(n) + i/D (f(n+1) - fiIm)] N D+) = 
=>} nD? f(n) + D(D-1)/2 > f(n) 
+ > n D(D-1)/2 (f(n+1)-f(n)) ] + (2D+1)(D+1)/6 Y 
(f(n+1)-f(n)) 
Come prima l'ultimo termine e' zero, e il secondo compensa Il 
terzo: 
=D*> nin) 


In conclusione MO = D }, f(n), e M®)= D2 5, n f(n). Attenzione 
che n rappresenta l'indice del tempo di campionamento in unita' 
di intervallo di campionamento. Ad esempio, con 
campionamento a 30 secondi, t = 2 h corrisponde a n = 2x120 = 
240. 


Immissione non istantanea 


Non e' strettamente necessario immettere il tracciante in un 
tempo molto breve. Supponiamo che l'immissione del tracciante 
sia descritta dalla funzione f(t), con f(t)=0 per t<0 e per t>t, e So” 
f(t) dt = 1, (normalizzata: f(t) ha dimensioni 1/S). 
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La curva di restituzione all'uscita risulta pertanto la convoluzione 
della risposta all'impulso h(t) con la curva di immissione, f(t), 
h'(t) = f f(t-s) h(s) ds 
dove l'integrale va da 0 a t. Quindi il momento di ordine zero 
Jh'(t) dt = f dt {-* ds f(t-s) h(s) 
=Sds h(s) [© dt f(t-s) 
=fds h(s) 
risulta invariato M', = Mg. Il momento di ordine uno risulta 
Sh'(t)tdt=fds h(s) f.® dt f(t-s) t 
=fSds h(s) So” dt f(t) (t + s) 
={dsh(s){s+t} 
dove t = f® t f(t) dt e' il momento di ordine uno della curva di 
immissione. Pertanto M', = M, + t M;. In conclusione con una 
immissione prolungata nel tempo Il momento di ordine zero non 
varia, mentre nella la formula del volume bisogna usare 
M'b9/(M'-TMo 


Per evitare che i residui non volatili non restino in grotta e' stato 
suggerito di spuzzare il tracciante su un canovaccio. Un 
semplice modello di questo modo di immissione e' un 
esponenziale decrescente, (1/t)exp(-t/t),. La formula della 
portata non varia, mentre nella la formula del volume bisogna 
usare 

M'b9/(M',-t Mg 


Un altro caso pratico consiste in un tracciante immesso in modo 
continuativo e costante a partire dal tempo t=0, cioe' f(t)=0 per 
t<0 e f(t)=k per t>0, come, per esempio, quando si ha una 
perdita di un inquinante. In questo caso h'(t) = fo! h(s) ds, 
pertanto h(t) = dh'(t)/dt. In questo caso h'(t) e' una funzione 
continua monotona non-decresente e per t molto grande tende 
al valore H = lim,_, h'0) = KSh(s) ds. 
Il momento di ordine zero di h(t) e' 

Mo = So” Ah'(8)/ds ds = h'(00) - h'(0) = H. 


1/ 


18 


Il momento di ordine uno risulta 
M, = Jo” Ss dh'(s)/ds ds = [s h'(S)]® - fo h'(8) ds 
= lim;_(TN'M-STN'(8) ds) 
= limp_- (So (H-h'(S)) ds) 
cioe' M, e' l'integrale della differenza fra H e h'(t). 


Con queste formule si possono quindi calcolare portata e 


volume a partire dalla curva di risposta. Da notare che il tempo 
iniziale di rilascio non appare in queste formule. 


Rami paralleli 


a g(z,t) 
Quando la grotta ha due ita 
rami paralleli parte del cn 
. > Vi, 2 
tracciante segue un agg” 


ramo e parte l'altro. Sia 
r(z,t) e'la z 
concentrazione del 
tracciante nal ramo di 
entrata, g(z,t) quella del primo ramo e (1-a) f(z,t) quella del 
secondo, dove a e' la frazione di tracciante che passa nel primo 
ramo e (1-a) nel secondo. 

Sg(z,t)dt = ff(z,t)dt = fr(z,t)dt 
Quindi attraverso una qualsiasi sezione trasversale della grotta 
(come il piano in figura) si ha un valore di MO costante: 

MO) = a fg(z,t)dt + (1-a) ff(z,t)dt = fr(z,t)dt. 
Pertanto la portata stimata da tutto il tracciante e' uguale a 
quella stimata dalla porzione di tracciante che passa in uno dei 
due rami: 

am/(afg(z,t)dt)=m/fr(z,t)dt. 


(1-a) f(z,t) \ 


Trasformata di Laplace 


Se si considera la grotta come formata da due (o piu') ambienti 
uno dopo l'altro, cioe' in serie, il momento di ordine zero e' lo 
stesso per entrambi gli ambienti (e per la grotta). Il momento di 
ordine uno della grotta e' la somma dei momenti di ordine uno 
dei due ambienti, cioe' e’ additivo. Anche p, e' additivo. 
L'input al primo ambiente e' un impulso (delta), ma l'input al 
secondo e' l'output del primo. Per ottenere la concentrazione in 
uscita al secondo ambiente bisogna fare la convoluzione delle 
risposte all'impulso dei due ambienti 

h(t) = Sh,(t-s) h,(s) ds 


| momenti possono esser espressi formalmente attraverso la 
trasformata di Laplace della funzione di risposta all'impulso, h(t), 
ovvero L[h(t)] = H(s) = {,” exp(-st) h(t) dt: 

Mo = H(0) 

M, = H'(0) / H(0) 

Ha = H"(0)/H(0) - {H'(0)/H(0)}£ 


Supponendo che la funzione di trasmissione di un tratto 
abbastanza uniforme di grotta dipenda dalla lunghezza, dalla 
velocita' di flusso, e dalla diffusione si puo' assumere che abbia 
due parametri: il ritardo di risposta, e l'ampiezza della risposta. 
Un percorso complesso puo' essere scomposto in tratti in 
parallelo e tratti in serie, come un circuito elettrico. Per due tratti 
in serie la funzione di risposta e' la convoluzione delle due 
funzioni di risposta 

f(t) = f,(t) * f.(t) =Jf,(t-s) f.(s) ds 
Per i tratti in parallelo, la somma 

f(t) = f,(t) + f.(t) 
Nei punti in cui un ramo si divide in due f, (t) + f.(t) = f(t) pertanto 
si puo' descrivere il bivio con un coefficiente di partizione a: f,(t) 
=af(t), e f.(t) = (1-a) f(1). 
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Passando alle trasformate di Laplace abbiamo che si 
moltiplicano per rami in serie e si sommano per rami in parallelo. 
Per esempio la grotta in figura ha risposta 

H(H+H)H, 


Una soluzione analitica 


Nell'ipotesi che la concentrazione dipende solo da z e t, 
l'equazione di diffusione diventa 

Oc/ot = - v Ac/oz + D dc/dz? 
dove si e' supposto anche che la velocita', v, e il coefficiente di 
diffusione, D, siano costanti. Questo modello si applica alla 
diffusione in un flusso a velocita' uniforme in un tubo. 
Le condizioni al contorno 

0c/or = 0 per r=R 

c(z,t) finito per r=0 
sono automaticamente soddisfatte poiche' abbiamo supposto 
che la concentrazione c sia funzione solo di z e t. 
La condizione iniziale e' un impulso al tempo 0 localizzato sul 
piano z=0 

c(0,t) = ò(t) 


Una soluzione e' 
c(z,t) = At!2 exp( - (z -vt)? / (4Dt) ) 
Infatti 
Oc/ot = { - 1/(2t) + v (z-vt)/(2Dt) + (z-vt)?/(4Dt?) } c 
v dc/oz = { - v (z-vt)/(2Dt) } c 
D dc/az? ={-1(20 + (z-vt)?/(4Dt?) } c 


Per verificare la condizione iniziale occorre calcolare 
limo (1/Vt) exp( - (Z -vt)? / (4Dt) ) 
Sez>QOsiha 


pal 


lim,_o €Xp(-z7/4Dt)/Vt = lims — co Vs exp( -s z2/4D ) = 0 
Sez=0siha 

lim,_o EXP(-v2t/4D)/Vt = 00 
Si calcola A in modo che la condizione iniziale sia soddisfatta, 

AS, lim,_o EXp(-‘2t/4D)/Vt dt = 2 A V(4D)/v fo? exp(-57) ds 
quindi 

A=v/V(4TD). 


La figura sotto (generata con wolfram-alpha) mostra come la 
concentrazione varia nello spazio e nel tempo. La 
concentrazione si sposta in avanti col tempo, e al contempo si 
allarga. 
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